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Observing dark energy dynamics is the most important aspect of the current dark energy research. 
In this paper we perform a global analysis of the constraints on the property of dark energy from 
the current astronomical observations. We pay particular attention to the effects of dark energy 
perturbations. Using the data from SNIa (157 "gold" sample), WMAP and SDSS we find that the 
best fitting dark energy model is given by the dynamical model with the equation of state across -1. 
Nevertheless the standard ACDM models are still a very good fit to the current data and evidence 
for dynamics is currently not very strong. We also consider the constraints with the recent released 
SNIa data from SNLS. 



I. INTRODUCTION 

In 1998 the analysis of the redshift-distance relation of 
Type la supernova (SNIa) has estabhshed that our uni- 
verse is currently accelerating P, ■ Recent observations 
of SNIa have confirmed the accelerated expansion at high 
confidence level0,0,|3. The nature of dark energy(DE), 
the mysterious power to drive the expansion, is among 
the biggest problems in modern physics and has been 
studied widely. A cosmological constant, the simplest 
candidate of DE where the equation of state (EOS) w re- 
mains -1, suffers from the well-known fine-tuning and co- 
incidence problems 0, Q- Alternatively, dynamical dark 
energy models with the rolling scalar fields have been 
proposed, such as quintessence J, 9], the ghost field of 
phantom [l^l and the model of k-essence which has non- 
canonical kinetic term [Til IT^. 

Given that currently we know very little on the theoret- 
ical aspects of dark energy, the cosmological observations 
play a crucial role in our understanding of dark energy. 
The model of phantom has been proposed in history due 
to the mild preference for a constant EOS smaller than 
— 1 by the observations [loj . Although in this scenario 
dark energy violates the weak energy condition(WEC) 
and leads to the problem of quantum instabilities 13], it 
remains possible in the description of the nature of dark 
energy [lj|. 

An intriguing aspect in the study of dark energy is 
that the recent SNIa observations from the HST/GOODS 
program^, combined together with the previous super- 
nova data, somewhat favor the dynamical dark energy 
model with an equation of state gett ing across -1 dur- 
ing the evolution of the universe 

dliilillli Ai- 

though the conventional scalar dark energy models also 
show dynamical behaviors with redshift, due to the in- 
stabilities of perturbations they cannot preserve the re- 
quired behavior of crossing the cosmological constant 
boundary ^19, .20, 21j. The required model of dark en- 
ergy has been called as quintom|0| in the sense that its 
behavior resemble the combined behavior of quintessence 



and phantom. However, the quintom models can be very 
different from the quintessence or phantom in the deter- 
mination of the evolution and fate of the universe ,2j|. 
There are a lot of interests in the literature recently in 
the building of quintom-like models. With minimally 
coupled to gravity a simple realization of quintom sce- 
nario is a model with the double fields of quintessence and 
phantomO HI Hi HI HI IH El . In such cases quin- 
tom would typically encounter the problem of quantum 
instability inherited in the phantom component. How- 
ever in the case of the single field scalar model of quin- 
tom, Ref. added a high derivative term to the kinetic 
energy and its energy-momentum tensor is equivalent to 
the two-field quintom model. Such a model with a high 
derivative term is possibly without quantum instabilities 
and as indicated by the SNIa observations, we are living 
with the ghosts 29]^. Perturbations of the quintom-like 
models have been studied extensively in Ref.|2C|. 

Observing the dark energy dynamics is currently 
the most important aspect of the dark energy study. 
Besides the SNIa data, a thorough investigation de- 
mands a fully consistent analysis of Cosmic Microwave 
Background(CMB), large scale structure(LSS) and with 
multi-parameter freedoms. The aim of current paper is 
to study the observational implications on dark energy 
in the consistent way. We extend our previous work of 
Ref. "2^ and study the full observational constraints on 
dynamical dark energy. In particular we pay great at- 
tention to the effects of dark energy perturbations when 
the equation of state gets across -1. Our paper is struc- 
tured as follows: in Section II we describe the method 
and the data; in Section III we present our results on the 
determination of cosmological parameters with the first 
year Wilkinson Microwave Anisotropy Probe (WMAP) 
pi, SNIa 4, 33] and Sloan Digital Sky Survey (SDSS) 
|34| data by global fittings using the Markov chain Monte 



^ For another interesting single field quintom model building see 
|3l. 



2 



Carlo (MCMC) techniques |35|, 133,133; finally we present 
our conclusions in Section IV. 



a divergence. In our study for such a regime, we match 
the perturbation in region 2) to the regions 1) and 3) at 
the boundary and set|2(l| 



II. METHOD AND DATA 

In this section we firstly present the general formula 
of the dark energy perturbations in the full parameter 
space of w especially when it crosses over the cosmologi- 
cal constant boundary. In our MCMC fittings to WMAP, 
SNIa and SDSS, we adopt a specific parametrization of 
the equation of state. 

Despite our ignorance of the nature of dark energy, it is 
natural to consider the DE fluctuation whether DE is re- 
garded as scalar field or fluid. In the extant cases like the 
two-field-quintom model as well as the single field case 
with a high derivative term (20.] , the perturbation of DE is 
shown to be continuous when the EOS gets across -1. For 
the conventional parameterized equation of state one can 
easily reconstruct the potential of the scalar dark energy 
if the EOS does not get across -1. Resembling the multi- 
field model of quintom or quintom with high derivative 
terms, the potential of the quintom dark energy can be 
directly reconstructed from the parameterized EOS on 
either side of the cosmological constant boundary^. In 
this paper we give a self-consistent method to handle the 
perturbation in all the allowed range of EOS especially 
for the region where EOS evolves close to and crosses -1. 

For the parametrization of the EOS which gets across - 
1, firstly we introduce a small positive constant e to divide 
the full range of the allowed value of the EOS w into 
three parts: 1) w > — 1 + e; 2) — 1 -I- e > w > — 1 — e; 
and 3) w < — 1 — e. Working in the conformal Newtonian 
gauge, one can describe the perturbations of dark energy 
as follows 381: 



S^O , 61 = 0. 



(3) 
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Neglecting the entropy perturbation contributions, for 
the regions 1) and 3) the EOS is always greater than —1 
and less than —1 respectively and perturbation is well 
defined by solving Eas. (|ll2|l . For the case 2), the pertur- 
bation of energy density S and divergence of velocity, 9, 
and the derivatives of S and 9 are finite and continuous for 
the realistic quintom dark energy models. However for 
the perturbations of the parametrizations there is clearly 



^ Although the multi-field dark energy models are more challeng- 
ing on theoretical aspects of naturalness, given that we know very 
little on the nature of dark energy, the energy momentum of such 
models can be identified wit h sin gle field scalar dark energy with 
high derivative kinetic terms Isflll . Our phenomenological formula 
of perturbations on DE corresponds to such models of multi-field 
(qui ntom ) with a negligible difference around the crossing point 
of -1I2I. 



In our numerical calculations we've limited the range to 
be I Aw = e| < 10~^ and we find our method is a very 
good approximation to the multi-field quintom, with the 
accuracy being greater than 99.999%. 

In the study of this pap er the parameterized EOS of 
dark energy is taken bv[39j| 



w{z) ~ Wq + Wi 



l + z 



(4) 



The method we adopt is based on the publicly available 
Markov Chain Monte Carlo package cosmomc^3i l4l|| . 
which has been modified to allow for the inclusion of dark 
energy perturbations with EOS getting across -l|i20j. We 
sample the following 8 dimensional set of cosmological 
parameters: 

P EE {iVb, i^c,Qs,T, Wo, Wl, Us, \og[lO'^° As]) (5) 

where ujh — flhh^ and uJc — flch^ are the physical baryon 
and cold dark matter densities relative to critical density, 
Qs is the ratio (multiplied by 100) of the sound horizon 
and angular diameter distance, r is the optical depth, 
is defined as the amplitude of initial power spectrum and 
rig measures the spectral index. Basing on the Bayesian 
analysis, we vary the above 8 parameters fitting to the 
observational data with the MCMC method. Through- 
out, we assume a fiat universe and take the weak priors 
as: r < 0.8, 0.5 < < 1.5, -3 < wq < 3, -5 < wi < 5, 
a cosmic age tophat prior as 10 Gyr< tg <20 Gyr. Fur- 
thermore, we make use of the HST measurement of the 
Hubble parameter Hq = lOO/i km s"^Mpc"^ 1^ by 
multiplying the likelihood by a Gaussian likelihood func- 
tion centered around h — 0.72 and with a standard de- 
viation a = 0.08. We impose a weak Gaussian prior on 
the baryon and density fJfc/i^ = 0.022 ± 0.002 (1 a) from 
Big Bang nucleosynthesis |43|. 

In our calculations we have taken the total likelihood 
to be the products of the separate likelihoods of CMB, 
SNIa and LSS. Alternatively defining = — 21og£, we 
get 



\total 



2 

XCMB 



XsNIa 



Xlss 



(6) 



In the computation of CMB we have included the first- 
year temperature and polarization data |3^ l44j with 
the routine for computing the likelihood supplied by the 
WMAP team . In the computation of nonlinear evo- 
lution of the matter power spectra we have used the code 
of HALOFIT [ill and fitted to the 3D power spectrum of 
galaxies from the SDSSfl^l using the code developed in 
Ref. III the calculation of the likelihood from SNIa 

we have marginalized over the nuisance parameter. For 
the main results of the current paper the supernova data 
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we use are the "gold" set of 157 SNIa published by Riess 
et aZ in In addition, we also consider the constraints 
from the distance measurements of the 71 high redshift 
type la supernova discovered during the first year of the 
5-year Supernova Legacy Survey (SNLS)j33j|. 



III. RESULTS 

In this section we present our results, particularly fo- 
cusing on the effects of the dark energy perturbation, we 
start with the descriptions on the background parame- 
ters, then present the constraints on dark energy param- 
eters. At last we give our constraints on dark energy 
from the recent observational data of SNLS. 

In Table 1 we list the mean la constraints on the pa- 
rameters with and without DE perturbations. We find 
that almost all the cosmological parameters are well de- 
termined in our both cases. However the reionization 
depth seems to be an exception where a vanishing r can- 
not be ruled out. We notice the determination on r is 
prior dependent. For example the WMAP collaboration 
have taken a prior like r < 0.3^^ lisl |4^ . which leads 
to a relatively stringent constraint on t by the observa- 
tions and a nonzero r is particularly favored by the high 
power of temperature-polarization cross correlation on 
the largest scales 50]. However when the strong prior on 
T is dropped one will in general get a less stringent bound 
from the full observational constraints, as also shown in 
Ref . ^3 • The prior on r is somewhat crucial for our pa- 
rameter estimation since its effects on CMB can be com- 
pensated with the tilt of the primordial scalar as well as 
the tensor spectrum. As we will show below it will also 
be correlated with the dark energy parameters due to the 
Integrated Sachs- Wolfe (IS W) eflPects. 

In Fig^ also we delineate the corresponding posterior 
one dimensional marginalized distributions of the cos- 
mological parameters from the combined observations of 
WMAP, SDSS and SNIa. The dotted vertical lines shows 
the quantity of every parameter with (red) and with- 
out (blue) DE perturbation giving the maximum likeli- 
hood. Due to the fact that the peaks in the likelihood are 
different from the corresponding expectation values, the 
dashed lines in Fig^ do not lie at the center of the pro- 
jected likelihoods. A vanishing r cannot be excluded by 
the full combined observations at high confidence level. 
For the order of parameters listed in (5) the best fit val- 
ues constrained by the full dataset (WMAP -I- SNIa -I- 
SDSS) is p=(0.023,0. 12, 1.04,0. 16,-1. 30, 1.25,0.995,3. 23). 
And for comparison the resulting parameters when 
switching off dark energy perturbations are given by 
p=(0.023,0. 12, 1.05,0. 14,-1. 15,0.63,0.962,3. 14). Compar- 
ing with the bottom of Table 1, although the minimum 
values have not been affected significantly (up to 1) 
by dark energy perturbations, all the best fit parame- 
ters have been changed. Moreover, the allowed param- 
eter space has been changed a lot and the constraints 
on the background parameters have been less stringent 



when including the dark energy perturbations. This can 
also be clearly seen from the two dimensional contour 
plots on the background parameters in FigI21 The rea- 
son is not difficult to explain. The ISW effects of the 
dynamical dark energy boosts the large scale power spec- 
trum of CMBpoj|. For a constant equation of state Ref. 
|5ll | has shown that when the perturbations of dark en- 
ergy have been neglected incorrectly, a suppressed ISW 
will be resulted for quintessence-like dark energy and on 
the contrary, an enhanced ISW is led to by phantom- 
like dark energy. In this sense if we neglect dark energy 
contributions, there will be less degeneracy in the deter- 
mination of dark energy as well as the relative cosmolog- 
ical parameters. However, dark energy perturbations are 
anti-correlated with the source of matter perturbations 
and this will lead to a compensation on the ISW effects, 
which result in a large parameter degeneracv|5lj. In fact 
as we have shown that crossing over the cosmological con- 
stant boundary would not lead to distinctive effects 2^ , 
hence the effects of our smooth parametrization of EOS 
on CMB can also be somewhat identified with a constant 
effective equation of state 52| 

_ Jdan{a)w{a) 

however the SNIa and LSS observations will break such a 
degeneracy. Thus for the realistic cases of including dark 
energy perturbations, the correlations between the dark 
energy and the background parameters as well as the auto 
correlations of the background cosmological parameters 
have been enlarged, as can be seen from Fig|21 

The contribution of dark energy perturbation affect 
significantly the distribution of wq and wi, which can 
also be seen from Fig|51on the constrains in the {wq, wi) 
plane. For the parameters {wq, wi) the inclusion of the 
dark energy perturbation change its best fit values from 
(-1.15, 0.63) to (-1.30, 1.25). In FigEl from out- 
side in, the contours shrink with adding 157 SNIa data 
provided by Riess et al and SDSS information. Dark en- 
ergy perturbation introduces more degeneracy between 
Wq and wi thus enlarges the contours significantly. For 
the discussion on dynamical dark energy we have sepa- 
rated the space of wq — wi into four areas by the lines 
of Wo = —1 and wq + wi = —1. The areas will then 
represent the quintessence where the EOS is always no 
less than -1 (the area with w > 1 can be reached by 
quintessence with a negative potential|5^), Quintom A 
where w is phantom-like today but quintessence-like in 
the past, phantom where the EOS is always no larger 
than -1 and Quintom B where dark energy has w > —1 
today but < — 1 at higher redshifts. From the figure we 
can see that dynamical dark energy with the four types 
are all allowed by the current observations and Quintom 
A seems to cover the largest area in the 2-dimensional 
contours with all the data we used. 

As shown in Fig|31 that wq and wi are in strong cor- 
relations. The constraints on w{z) are perhaps rel- 
atively model independent, as suggested by Ref.0]. 
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TABLE 1. Mean la constrains on cosmological parameters using different combination of WMAP, SNIa and SDSS 
information with/witliout DE perturbation. For the weakly constrained parameters we quote the 95% upper limit instead. 





With DE Perturbation 


Witliout DE Perturbation 




WMAP 


WMAP±SN 


WMAP±SN±SDSS 


WMAP 


WMAP±SN 


WMAP+SN±SDSS 






0234+° °°^° 

'-'•"^'5^-0.0011 


0.0232 ± 0.0010 


0.0235 ± 0.0013 


0.0232 ±0.0011 


0.0230 ± 0.0009 




0.124 ±0.016 


0.128 ±0.018 


0.123 ±0.010 


0.111 ±0.020 


0.119 ±0.018 


0.122 ±0.010 


Qs 


1.046 ± 0.006 


1.047 ± 0.006 


1.046 ± 0.005 


1.046 ± 0.006 


1.046 ± 0.006 


1.046 ± 0.005 


T 
Wo 
Wi 
Us 

log[10'°4s] 


< 0.256(95%) < 0.264(95%) 

r, 700+0.623 -i 1 70+0-231 
"•'■J^-0.613 J^-J^ '^-0.226 

< 1.59(95%) 0.361115:^83 

077+0.029 Q QOC+0.030 
"•^''-0.030 ^-y^-COSl 

3.181 ±0.134 3.207 ±0.132 


< 0.256(95%) 

1-10'-0.190 
c:Q7+0-657 
'-'•^^'-0.713 

0.982 ± 0.030 


< 0.399(95%) 

-0 R17+0193 
U-Dl' -0.190 

< 0.832(95%) 

0.995+nS 
3.2451°;?°^ 


< 0.324(95%) 

-i.080l°:Jf7 

0.359l°;f7^ 
0-98ll°;gg 
3.206l°:l^g 


< 0.246(95%) 

1 008+0-078 

-i.uy»_o.o8o 
416+°-^^^ 

"•^l"-0.153 

0.97010;E;24 

+0.118 
' —0.119 




r, 7(.o+0.073 


0.678 ± 0.045 


0.68ltH^3i 


0.681 ± 0.086 


0.697 ± 0.048 


0.686 ± 0.031 


Age/GYr 


13.45 ±0.27 


1 ry+O.SO 


13.67 ±0.24 


13.60 ±0.30 


13.62 ±0.29 


13.67 ±0.23 




r, 207+0.072 
u.^itl _o.073 


0.322 ± 0.045 


n Q1Q+0 030 


0.319 ±0.086 


0.303 ± 0.048 


0.314 ±0.031 


Ho 


n 007+0.152 

^•y^'-0.154 
14 ■:{t;+4.71 
l*.O0_4.68 
71 71+7.95 


o.9i3ts;:ifo 

14.72tr* 


0.8541^^7 
14.4ltt^^ 
67.90l^:fg 


0.818 ±0.120 

17 1 7+6.94 


15.49t^:f3 
68.89l^;52 


n a»o+0-073 

13.73tr,« 
68.07t^;i 


xVd.o.f 


1429.1/1342 


1610.4/1499 


1633.2/1518 


1428.4/1342 


1610.6/1499 


1634.1/1518 



Following we obtain the constraints on w{z) by com- 
puting the median and 1, 2a intervals at any redshift. In 
Fig0|we plot the behavior of the dark energy EOS as a 
function of redshift z, we find that at redshift z = 0.3 
the constraint on the EOS is relatively the most strin- 
gent. One can see that the perturbation reinforces the 
trend of DE to cross -1 at z ~ 0.3. However due to the 
limitation of the observational data, the quintom scenario 
is only favored at Icr by the full dataset of WMAP, SDSS 
and SNIa. We find the value at z — 0.3 is restricted at 

= 0.3) = -i.oo2iEJ:°^ro:l|° (s) 

for the case without dark energy perturbations and 

w{z ^ 0.3) = -im9t"of,ltr2ll (9) 



when including dark energy perturbations. Correspond- 
ingly at redshift z = 1 the constraints turn out to be 



w{z 



1) 



-0.890 



+0.180+0.193 
0.159-0.301 



without perturbations and 

w{z = 1) = -0.868 



-0.215+0.520 
-0.204-0.815 



(10) 



(11) 



when including dark energy perturbations. One should 
bear in mind that such a constraint is not really model 
independent, as shown in Refs. '56^, (&3|- 

Recently the authors of Ref.,2^] made the distance 
measurements to 71 high redshift type la supernovae 
discovered during the first year of the 5-year Super- 
nova Legacy Survey (SNLS). SNLS will hopefully dis- 
cover around 700 type la supernovae, which is an in- 
triguing ongoing project. Following Ref.|3^ we combine 



the "new 71 high redshift SNIa data the 44 nearby 
SNIa", together with WMAP and SDSS. We plot the 
constrains on the dark energy parameters in Fig[Sl We 
find the current data of SNLS are very weak in the de- 
termination of dark energy parameters. Although our 
best fit values are given with a quintom like dark energy: 
{wq, Wi) — (—1.47,1.44), a cosmological constant fits 
well with SNLS in the Icr region.^ 




FIG. 1: 1-D constrains on individual parameters us- 
ing WMAP+157 "gold" SNIa+SDSS with our 8-parameter 
parametrization discussed in the text. Solid curves illustrate 
the marginalized distribution of each parameter with(red) and 
without(blue) DE perturbation. Dotted vertical lines shows 
the quantity of every parameter with (red) and without (blue) 
DE perturbation giving the maximum likelihood. 



3 Ref. 1581 made some study on SNLS implications of dynamical 
dark energy using SNIa data only. 
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FIG. 2: Sampled 2-D contours of the background parameters 
and also the contours among dark energy and the background 
parameters. Here we use the same parametrization and data 
sets as FIG. 2, red solid and blue dashed lines are for pertur- 
bated and unperturbated DE respectively. 




-2 -1 -2-10 



FIG. 3: 95% constrains in the {wo, wi) plane with and without 
dark energy perturbation from left to right. Shaded dark 
orange region is excluded by WMAP only for our 8-parameter 
estimation. The dashed lines stand for wo = —1 and wo + 
Wi = —1, see the text for details. 



IV. DISCUSSION AND CONCLUSION 
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FIG. 4: Constrains on w(z) using WMAP + 157 "gold" SNIa 
data -|- SDSS with/without DE perturbation. Median (central 
hne), 68%(inner, dark grey) and 95%(outer, light grey) inter- 
vals of w(z) using 2 parameter expansion of the EOS in (4). 
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™0 

FIG. 5: Two dimensional constraints on the parameters of dy- 
namical dark energy from the combined constraints of SNLS, 
WMAP and SDSS. The solid and dotted lines are the con- 
straints with and without dark energy perturbations, respec- 
tively. 



In this paper we have performed an analysis of global 
fitting on the EOS of dark energy from the current 
data of SNIa, WMAP and SDSS. Dark energy pertur- 
bation leaves imprints on CMB through ISW effects and 
changes the matter power spectrum by modifying the lin- 
ear growth factor as well as the transfer function. Our 
results show that when we include the perturbations of 
dark energy, the current observations allow for a large 
variation in the EOS of dark energy with respect to red- 
shift''. A dynamical dark energy with the EOS getting 



across —1 is favored at la with the combined constraints 
from WMAP, SDSS and the "gold" dataset of SNIa by 
Riess et al. When we use the recently released SNLS 
data and the nearby data of type la supernova instead, 
the parameter space is enlarged and a cosmological con- 
stant is well within Icr, although a quintom dynamical 
dark energy is still mildly favored. 

In our all results, the perturbation of dark energy plays 
a significant role in the determination of cosmological pa- 
rameters. Neglecting the contributions of dark energy 
perturbation will lead to biased results which arc more 



* Ref. l59l studied the perturbations of dy nam ical dark energy 
only for the regime where to > — 1 and Ref. l6r| considered both 
the cases for w > —I and w < —I, but did not include the per- 



turbation s for qu intom like dark energy; previous global analysis 
like Ref. I54ll6ll did not consider dark energy perturbations. 
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stringent than the real cases. In the next decade, there 
wiU be many ongoing projects in the precise determina- 
tion of cosmological parameters. We can hopefully de- 
tect the signatures of dynamical dark energy like quin- 
tom through global fittings to the observations, where 
it is crucial for us to include the contributions of dark 
energy perturbations. 
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